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Electron and phonon correlations in systems of one-dimensional electrons coupled to phonons are 
studied at low temperatures by emphasizing on the effect of electron-phonon backward scattering. 
It is found that the 2fcF-wave components of the electron density and phonon displacement field 
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coupled to one-dimensional or three-dimensional phonons, and they are long-ranged for repulsive 
electron-electron interactions for a three-dimensional array of parallel one-dimensional conducting 
chains coupled to three-dimensional phonons. 
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I. INTRODUCTION 

The physics of one-dimensional (ID) interacting elec- 
tron systems has attracted lots of interests. The low- 
energy long-wavelength excitations of such systems can 
be described by a Tomonnaga-Luttinger liquids (TLL) 
model [lj|. Such systems exhibit a critical like behav- 
ior (power law) of the charge-density correlations at zero 
temperature. Thus, ID electron systems are always at 
the verge of an instability without being to order [§] . On 
the other hand, it was noticed by Peierls Q that due 
to electron-phonon interactions a charge density wave 
(CDW) accompanied by a periodic lattice distortion ap- 
pears in the ground state of a ID metal, both periods 
being n/kp. For an arbitrary band filling the period 
is incommensurate with the underlying lattice, where 
the number of electrons per site is not a simple frac- 
tional number (eg., 1/2, 1/3, etc.). Later studies Q 
on incommensurate systems of ID electrons coupled to 
ID phonons (1D1D) show that the tendency to form 
CDW is increased significantly by electron-phonon cou- 
pling; however, true long-ranged-order CDW is still miss- 
ing. In reality ID electron systems are often embedded 
in three-dimensional (3D) environment, and phonons in 
such systems are effectively 3D. Therefore, it is interest- 
ing to know if the coupling between ID electrons and 3D 
phonons could lead to long-range-order CDW. 

For simplicity, we limit this paper to spinless incom- 
mensurate systems. Our work focuses on the effect of 
the electron-phonon coupling at large momentums (i.e., 
Backward-scattering process). We show that the Ikp- 
wave components of the electron density and phonon 
displacement field share the same correlations at large 
distance. This result applies to all the systems we con- 
sider in this paper. It is exact and does not depend on 
the electron-phonon coupling strength, phonon frequency 
and electron-electron interactions. This result implies 
that phonon correlations are quasi-long-ranged even for 
arbitrary weak electron-phonon coupling in 1D1D, in 
which previous studies Q showed that electron corre- 
lations are quasi-long-ranged. 

The first system we consider is a single conducting 



chain coupled to 3D acoustic phonons (1D3D). We find 
that electron-phonon coupling in this system does not 
change the qualitative power-law electron correlations, 
although it does increase the tendency to form CDW sig- 
nificantly (i.e., make electron correlations decay much 
slower). Surprisingly, our calculation shows that 1D3D 
even has less increased tendency to form CDW than 
ID ID, given that both the phonon frequency and the 
electron-phonon coupling strength are the same in these 
two cases. Furthermore, phonon correlations are also 
quasi-long-ranged in 1D3D. 

The second system we discuss is a 3D array of paral- 
lel ID conducting chains couple to 3D acoustic phonons 
(3D3D). For arbitrary weak electron-phonon coupling, 
the system undergoes a quantum phase transition driven 
by varying electron-electron interactions: for repul- 
sive electron-electron interactions, electrons on different 
chains are effectively correlated, which leads to long- 
ranged-ordered CDW and hence long-ranged phonon cor- 
relations; for attractive electron-electron interactions, 
electron correlations remain ID, and the system is effec- 
tively a simple sum of 1D1D, which leads to quasi-long- 
ranged electron and phonon correlations. In the limit 
of strong coupling and/or low phonon frequency, both 
electron and phonon correlations are long-ranged for re- 
pulsive electron-electron interactions, and interestingly, 
the phase fluctuations of the 2&F-wave components of 
the electron density and phonon displacement field are 
not independent but locked together. The 3D3D system 
has also been considered recently by Artemenko et al [8( . 
Their work agrees with ours in the limit of strong cou- 
pling and/or low phonon frequency. 

It was argued [4| that the electron-phonon coupling at 
small momentums (forward-scattering process) is not im- 
portant since its effect is of the order c 2 jv F <C 1 (vf is 
the Fermi velocity, and c is the sound velocity). How- 
ever, it was pointed out [5| that in strongly correlated 
systems coupling to small-momentum acoustic phonons 
can lead to the Wentzel-Bardeen (WB) singularity [f| @] 
and hence become important. This singularity is a crit- 
ical point at which for a critical electron-phonon cou- 
pling constant, the attractive electron-phonon interac- 
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tions cause the system to be unstable. It was shown [f| 
that the critical point can be reached as one approaches 
the half- filling for the Hubbard model [15j |. In this pa- 
per we assume that the system is far away from the WB 
singularity for a general filling, and the electron-phonon 
coupling at small momentums is thus not important. In 
fact we even find that the WB singularity does not exist 
in 1D3D no matter the filling. This is because the bulk 
phonon freedom suppresses phonon fluctuations on the 
chain, thus making the system robust. 

This paper is arranged as follows. In Sec. [IT] we intro- 
duce the model. In Sec lIIII we review 1D1D systems. In 
Sec. ITVland Sec. |V] we discuss 1D3D and 3D3D systems, 
respectively. Finally some technical details are given in 
the appendixes. 



II. MODEL 

In this section we discuss the model for 3D3D, from 
which the models for 1D1D and 1D3D can be obtained 
as special cases. Let us start with electrons. Using the 
TLL description for electrons 0, Q , the electron density 
inside the jth chain can be written in the form 
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cos (2k F + %<t>j 



(1) 



where <pj is a slowly varying phase, x denotes the direc- 
tion along the chains, a is the cut-off and of the order of 
the lattice spacing, kp is the Fermi momentum. The first 
term on the right-hand side of |T]) is the long-wavelength 
part of the electron density, the second term is the fast 
oscillating part. The total electron action is a sum of 
TLL models 
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where v is renormalized Fermi velocity, K is a dimen- 
sionless parameter, which is bigger than 1 for attractive 
electron-electron interactions and less than 1 for repulsive 
electron-electron interactions, r is the imaginary time. 

Since electrons are restricted on the chains, an effective 
phonon action in terms of the on-chain freedom is needed. 
It can be obtained by starting with the standard lattice 
version of a 3D phonon action and eliminating the off- 
chain phonon freedom. This effective action is expected 
to have the following form 

s Wi h = dk?/ dxdT H^) 2 

+K x a 2 {d x u j f + K ± \ Uj - u 3+l \ 2 } , (3) 

where h is Planck Constant, a is the lattice spacing, b 
is the distance between nearest-neighbor chains, m is 
the atom mass, _L denotes directions perpendicular to 
the chains, K x and K± are the effective spring constants 
along x and _L directions, respectively. 



In general the action for electron-phonon coupling can 
be expressed in the form [l(| 



Sep/H = ^ ^ J dxdr (d x Uj)pi(x) , 



(4) 



where 7 is the coupling constant. Plugging (TT]) into the 
above action we obtain 



Sep/K = J dxdr -(d x Uj)(d x <pj 

j 

cos (2k F x + 2cj)j)] . 
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The first piece on the right-hand side of ([5]) corresponds 
to the electron-phonon coupling at small momentums, 
which is responsible for the WB singularity. We will dis- 
cuss it in section [TV] The second piece is the coupling 
at large momentums, which is the primary focus of this 
paper. For this latter type of coupling, only the phonon 
field with wavelength near l/2kp is important due to the 
fast oscillation of the cosine term. We write 



1, 



u j (x) = ~[e 2ik ^ J (x)+H.c], 



(6) 



where ipi is a slow varying complex field. Then the action 
for the coupling at large momentums can be rewritten as 



S l ep /h — E J dxdT (^ikpipj + d x "ijjj 
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where the piece which involves e AlkF ' x vanishes when inte- 
grated over x, also d x ipj is negligible compared to 2kptpj 
at long wavelength (i.e., q x < 2kp). The action can be 
further simplified as 



(8) 



At this point, it is also convenient to have a phonon ac- 
tion in terms of ipj . Plugging ^ into ([3]) we obtain 

SfflhVh = 5^E/ dxdr{^\d T ^\ 2 + 2K x a 2 k F x 



\^\ 2 + ^ 2 \dM 



i+H 



(9) 



For convenience, we rescale the lengths and field such 
that the actions are fully expressed in terms of dimension- 
less quantities. An appropriate rescaling is the following 



cit/(ttv), x—tax/n, 
brj_/ir 



(10) 
(11) 
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After the rescaling, the actions ((2]), ([9]) and ([8]) become 
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where Cj_ = K±a 2 /mv 2 , the filling factor n e = akp/ir 
(i.e., the averaged number of electrons per lattice site) 
which is of order 1 for a general filling. Both g\ and gi 
are dimensionless quantities defined as 
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where c = a^J K x /m is the sound velocity along the 
chains, 0J2k F is the phonon frequency at wave vector 2kp, 
h 2 /ma 2 is approximately the ground state energy of an 
atom with mass m trapped in an infinite one-dimensional 
potential well with width a, Hirv/a is of the order of the 
Fermi energy. Typically g\ is of order 10 -5 — 10~ 6 . In 
the present paper we assume 
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which can be fulfilled for a general filling. This condi- 
tion actually implies that the system is far away from 
the WB instability, so that the electron-phonon coupling 
at small momentums is not important. For the Hub- 
bard model very close to half-filling, since g\ /q 2 which is 
proportional to v 1 ^ 2 drops rapidly to zero 15 1, the pos- 
tulates (fT"7)) breaks down; therefore, the electron-phonon 
coupling at small momentums could be important. This 
will be further discussed in section HVl 



III. 



A SINGLE CONDUCTING ID CHAIN 
COUPLED TO ID PHONONS 



The model for 1D1D can be recovered by restricting 
the index j to be 1 in the actions (fl"2"j) . (fl"3"j) . fT4")) and 
setting C± to be 0. To lighten the notation we drop the 
index j in this section. First we want to show a very 
general result which is exact and does not depend on the 
values of K, the electron-phonon coupling strength and 



phonon frequency. We introduce a new complex field , 
which, in Fourier space, is related to </> and ip by 

ip' R (q) = MO) ~ 292 °^ [ dxdr sin 20 e^, 



(18) 

i>'M = Mo) - 292 °^ j dxdr cos 20 



(19) 

where V = Lhv/{k B Ta 2 ), G^(q) = q 2 T + g 2 (l/2n e ) 2 q 2 x + 
g 2 , the subscripts R and / denote the real and imaginary 
components, respectively. In terms of and 4> the total 
action can be nicely separated into two parts which are 
completely decoupled from each other: 



(20) 
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[ dxdrdr' G vH (t - t') cos [24>(x,t)- 
gi J 

20(x,r')], (21) 
where the phonon correlation function is given by 

G ph (r) «e-^M. (22) 



To obtain (f2Tj) we have neglected the phonon dispersions 
(d x 4>) 2 since the leading order effect on electron correla- 
tions comes from the phonon kinetic energy (d T <f) 2 . Since 
ijj' and (f> are decoupled, from (fl8|) and (|19p the fluctua- 
tions of ifj can be calculated as 

^{q)i>{-q)) = IGt: 2 a 2 g 2 2 G 2 ph (q) J dxdr [e^x 

(PtoAVPtoAfyil+Gph®, (23) 

where pik F is the 2/c^-wave component of the electron 
density. Since G p h (q) tends toward a constant for q <C g±, 
this result implies that at large distance (i.e., r » g^ 1 ) 
the 2fc^-wave component of the phonon displacement 
field has the same correlations as pik F ■ 

To understand the detailed behavior of the correla- 
tions, more calculations are needed. In the limit of weak 
coupling and/or high phonon frequency (i.e., g 2 g^ 2 <C 
1), the effort is put on the phonon- mediated effective 
electron action given by (|2ip. Since the exponentially 
decaying phonon correlation G p h{j) imposes a cutoff 
at r ~ at the wavelength longer than l/<?i the 

model ([2Tj) essentially reduces to the one for free elec- 
trons but with renormalized K, which we denote as K' . 
A Gaussian variational method [131) and a perturbative 
renormalization group (RG)d 1 1 31 ] have been applied to 
this model. Both find that electron-phonon coupling in- 
creases the tendency to form CDW. Specifically, in terms 
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of K' the result is expressed as 



using a Gaussian variational method, we get 
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The details of the calculation are given in appendix [X] 
To obtain the above result, it has been assumed that the 
renormalization of K is small, that is, (K/K 1 ) 2 — 1 <C 1, 
which leads to a self-consistent condition. For K < 1, 
this self-consistent condition is just 



^«9\- K - 
9i 



(25) 



which prohibits the application of the Gaussian varia- 
tional result in the region above the locus OA in Fig. 
® 

When the condition (|25|) is violated, neither the Gaus- 
sian variational method nor the perturbative RG can give 
an analytical result. A new strategy is used. A mean- 
field theory shows that ip saturates for K < 1. There- 
fore, it is believed that for K < 1 and in the limit of 
strong coupling and/or low phonon frequency, the im- 
portant fluctuations of ip are its phase fluctuations, and 
the amplitude fluctuations are massive and hence negli- 
gible. Based on this belief we now work on the truncated 
action: 
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where ipo is fixed at the saturated value, and the con- 
stant pieces have been thrown away. Hereafter we name 
this strategy as "fixed-amplitude approximation" just for 
convenience. In the following we will first estimate the 
saturated value of ipo, then calculate the fluctuations of 
and 9 from the truncated model (I26|) . and finally discuss 
the limits of the validity of the fixed-amplitude approxi- 
mation. 

We estimate the saturated value of if> using a mean- 
field theory. Approximating ip as uniform both in space 
and time by doing = if}Q, the total action reduces to 

SlZjh =\J dxdr -L [{drcpf + [d^f] 
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where the electron part of the action is just a Sine-Gordon 
model. Directly integrating out the field in favor of ipQ 
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The detailed calculation is shown in appendix [Bl The 
first term in the action [)28p is negative, which means that 
a periodic lattice distortion with period 7t//cf lowers the 
electron energy; the second term is the elastic energy for 
such a lattice distortion and hence positive. To find the 
ground state, we minimize the action (|28j) with respect 
to ipQ. We find 



■00 = 9i 



(29) 



for K < 1 and ipo — for K > 1. This indicates that 
a periodic lattice distortion with period 7r/fci? appears 
in the ground state for repulsive electron-electron inter- 
actions; however, -0 is not long-range-ordered due to its 
phase fluctuations. Later in this section we will show 
that tp is quasi long-range-ordered with power-law cor- 
relations. For K not very close to 1, the result given 
in (|29[) qualitatively agrees with that in ref. [3J], in 
which the electron freedom is eliminated by mapping the 
model (|27p to an exactly solvable field-theoretical model 
(i.e., Massive Thirring model) and calculating the ground 
state energy of electrons. In the limit of K — > 1~, (j2"9")l 
shows 0>o - * 0, and the classic Peierls theory predicts 
■00 ~ exp (— gf/g%) /<?2, which is essentially very small. 

The above mean-field calculation assumes that the sin- 
gle electron gap opened at the Fermi level is much less 
than the Fermi energy, that is, {g 2 i>o) 1/{2 ~ K) < 1. This 
assumption combined with (f2"9"| leads to the postulated 
condition (fT7)) . or equivalently 



01*00 < 1 



(30) 



Now we discuss the truncated model (|2"6"|) . This model 
has been studied previously in the context of liquid crys- 
tal [Tl|. It can be shown that the coupling term in this 
model is either relevant or irrelevant depending on the 
values of ipo and K. The critical values of tpo and K can 
be easily calculated in the limit K — > (i.e., the electron 
field is locked) or 0o — > oo (i.e., the phonon phase field 
9 is locked) since in both limits the model reduces to the 
Sine-G ordon m odel. We get K c = 2 for ipo — > oo and 
iJjq — \Jn e /irgi for K — > 0. In general, the critical values 
of 0>o and K satisfies 
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The detailed derivation is given in appendix [C] This 
function defines a critical line in the iJ)q 2 —K plane, which 
is illustrated in Fig. [TJ Thus, for K < 1, only for which 
we can get a non-zero ipo, we have 
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FIG. 1: The line EF is denned by the function given in (|31[1 . 
The cosine term in the action (|26p is irrelevant above EF and 
relevant below it. 



For tpQ < tpQ, the cosine term becomes irrelevant, and <f> 
and 9 are decoupled. For ip$ > ipQ, the cosine term is rele- 
vant, and the fluctuations of cj> and 9 are bound together. 
This actually implies a self-consistent condition for the 
fixed-amplitude approximation; that is, the estimated ipo 
must be bigger than ?/>q such that phase fluctuations do 
not drive the cosine term into irrelevance. This cosine 
term offsets the positive phonon energy and is necessary 
for obtaining a non-zero tf>o m the first place. In terms 
of <7i and gi this self-consistent condition is given by 



9 2 ^ - 
— > Si 
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(33) 



which prohibits the using of the fixed-amplitude approx- 
imation in the region below the locus OB in Fig. [2] For 
the systems satisfying (f3"3"|) both the effective actions of (f> 
and 9 can be obtained by making 2</> = 9 in the truncated 
model (f!26"|) , and they are given by 

S[E]/H = J dxdr[R(d T Ef+B(d x E) 2 ] ,(34) 

where S can be either 4> or 9/2, and the parameters R 
and B are 



R = {l + AirKift), 



B = 1 + nKgl 
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(35) 
(36) 



This result is consistent with our earlier conclusion that 
at large distance the 2kp -wave components of the phonon 
displacement field and electron density share the same 
correlations. Note that the correction to B comes from 
the phonon dispersions along the chains and is much less 
than 1 according to the assumption (|30[) . Also it is of 
the order of the correction to R multiplied by g\ and 
hence negligible. Plugging (|2"5|) into ([33)1 and defining 



FIG. 2: Limits of the validity of the two approaches. The re- 
sult of the Gaussian variational method does not apply in the 
region above OA, while the fixed-amplitude approximation 
does not apply in the region below OB. The regions below 
OA and above OB are defined by the conditions (|25|l and 
([33)1 in the ID ID case and by (J7HJ and (86) in the 3D3D case, 
respectively. In the intermediate region between OA and OB 
neither of the two approaches applies. 



K' = K/VR, we get 
K> 



fJi 



(37) 



Now let us discuss the validity of the fixed-amplitude 
approximation. In addition to the self-consistent condi- 
tion J33|) caused by phase fluctuations, there is another 
one caused by amplitude fluctuations, which is 
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where <5"0o is the deviation of the amplitude of ip from the 
saturated value. For a crude estimation of ((Sipo) 2 ), we 
assume that the effective action for Sipo (to the quadratic 
order) has the following form 



S[5f }/h 
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where the mass term is obtained by expanding the mean- 
field action (|28[) around its minimum, and the derivative 
terms come from the free phonon action. Then after a 
straightforward calculation, the condition (l38|) also leads 
to (I33[) coincidentally. 

In both limits the charge-density correlations can be 
calculated, given by the non-universal power laws 



(P2k F {r)p2k F (0)) - r 



-IK 1 



(40) 



where r = (x,t). This implies that in electron-phonon 
coupling does not change the qualitative TLL behavior, 
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although the tendency to form CDW is increased since 
K' is smaller than K. Based on the general conclusion 
given in (|23j). at large distance phonon correlations are 
also governed by the same non- universal power laws: 



(^(f)V(O)) ~ r 



-IK 1 



(41) 



Surprisingly, this implies that in 1D1D while electron- 
phonon coupling does not have a qualitative effect on 
electron correlations, it does change phonon correla- 
tions qualitatively, that is, the correlations of the 2k p- 
wave component of the phonon displacement field, which 
are short-ranged in the absence of electron-phonon cou- 
pling, become quasi long-ranged. According to the self- 
consistent conditions given by ([25]) and ([33]) . there is 
an intermediate region in which neither of the two ap- 
proaches apply, and therefore no quantitative result is 
available. This region is illustrated in figure[2]as the area 
between the loci OA and OB. However, in this region 
we do expect the same qualitative conclusion. 



In the limit C± — > 0, the renormalized ID phonon prop- 
agator given by (|45|) reduces to the bare ID one (i.e., 
the one which is not renormalized by the bulk phonon 
freedom) as expected. For Cj_ being non-zero, for small 
qs the renormalized one has a logarithmic dependence 
on q, while the bare one has a second-order power-law 
dependence. This difference has a drastic effect on the 
electron-phonon coupling at small momentums. We de- 
lay this discussion to the later of this section. Now we 
focus our attention on the electron-phonon coupling at 
large momentums (i.e., near 2kp). Expanding the action 
(|4~4"]) around q x — ±2akp/7r, q T = (the factor a/ir is due 
to the rescaling), we get 



2akj. 



(46) 



where the complex phonon field iji is the slow varying part 
of the 2/c^-wave component of the displacement field u' 



IV. A SINGLE CONDUCTING ID CHAIN 
COUPLED TO 3D PHONONS 

In this section we consider a single ID conducting chain 
coupled to 3D phonons. In this case the action for elec- 
trons and electron-phonon coupling is the same as that 
in 1D1D, and the phonon action is given in ([3]). For con- 
venience, we transform the phonon action into its contin- 
uum version and perform the dimensional rescaling given 
in (fTUl [TT|) and the field rescaling 
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After this procedure, we obtain 



S$[u']/h = ±j d 3 rdr 
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Since electron-phonon coupling only happens on the 
chain, an effective phonon action for the on-chain free- 
dom is needed. We integrate out the bulk off-chain 
phonon freedom in favor of the on-chain one u' (x, r) = 
u'(x,t, fj_ = 0), where we have assumed that the posi- 
tion of the chain is at r± = 0. This leads to an effective 
ID phonon action 
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where 
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The inverse of the propagator in (|46| , to the lowest order 
of q x ,, is given by 



G" 



with 
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n 9193 
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9f + *C ± /n 



(48) 

(49) 
(50) 
.(51) 



To get (|48|) , we have thrown away the term linear in 
q x (i.e., (fel^l 2 )! which only leads to boundary effects. 
Plugging (gHJ) into (gBJ we get 



dxdr 



9z\d^f 



9A\d x ip\' 



(52) 



Therefore, with respect to the electron-phonon coupling 
at large momentums, the effective model of phonons in 
1D3D is qualitatively the same as that in 1D1D but with 
modified coefficients. Thoughtful reader might have seen 
this right from the beginning [12]. Note that g[, 53 and 54 
are all increasing functions of C± . This is expected since 
the bulk phonon freedom suppresses phonon fluctuations 
on the chain. 

Since we have shown that the models for 1D3D and 
1D1D are essentially the same, the calculations should be 
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also very similar. Repeating the (virtually same) calcula- 
tion, wc obtain, for 1D3D and repulsive electron-electron 
interactions, 

in the limit of weak coupling and/or high phonon fre- 
quency and 

-l~S3(»2»i J (54) 

in the limit of strong coupling and/or low phonon fre- 
quency. It can be verified that both K' obtained from the 
two limits are increasing functions of C±. Therefore, we 
conclude that a single ID chain coupled to 3D phonons 
does not lead to long-range-ordered CDW, instead the 
increased tendency to form CDW in such systems is even 
less than that in ID ID, given that the two dimension- 
less quantities g\ and gi (or equivalently, the phonon fre- 
quency and the coupling strength) are the same in these 
two cases. This result does not seem obvious just from 
qualitative arguments. On one hand, the bulk phonon 
freedom increases the energy cost to form a lattice dis- 
tortion along the chain, which is manifested in that g[ is 
a increasing function of C± . This is a negative effect on 
the increased tendency to form CDW. On the other hand, 
both <73 and g^ are increasing functions of C± , which en- 
hances the correlations of ip. This is a positive effect on 
the increased tendency to form CDW. Only a quantita- 
tive calculation can tell that the former effect dominates. 
In principle we could also think about a single ID chain 
coupled to 2D phonons. We expect that in such sys- 
tems long-range-ordered CDW is still missing, and the 
increased tendency to form CDW is less than that in 
1D1D but more than that in 1D3D, given that both g\ 
and (/2 are the same in the three cases. In addition, the 
2fcp-wave components of the phonon displacement field 
in 1D2D and 1D3D are quasi long-range-ordered due to 
electron-phonon coupling. 

Now we discuss the electron-phonon coupling at small 
momentums. According to the action given in (|5|), the 
action for the electron-phonon coupling at small momen- 
tums is given by 



S s ep /h= J dxdr (d x u )(d x 



(55) 



After the dimensional and field rescaling which are given 
in (flO|) and (|42f respectively, the above action becomes 

si /n = J^2l f dxdT (d x u' )(d x 4>) . (56) 



The effective action for the on-chain phonon field is given 
by (HU). Since the total action is quadratic in u' , we 
integrate it out. This leaves us an effective electron action 



ototal 
0„; 



2irK 



E 



8? 



^k(^) q 

\ an e I 



xcf>(-q)<j>(q) , 



(57) 



where G(q) is the phonon propagator. The first two terms 
on the right-hand side of ([57)) are from free electrons, and 
the third term, which is retarded and negative, is con- 
tributed by phonons. First let us treat phonons as pure 
ID. Plugging G(q) as the bare ID phonon propagator 
into (|57|) . we get 



5*. 



total 



i 



2izK 



E 



2^K^g\q x 



X<f>(-q)(f>(q) . 



(58) 



In the limit q T — > 0, the retarded part is proportional to 
q x and cancels the first term at the point 



1 

K 



8tt j 



(59) 



at which the electron density becomes unstable towards 
long wavelength fluctuations. This singular point is re- 
ferred to as the Wentzel-Bardeen (WB) singularity. Since 
K and a/ a are both of order 1, as long as the postulates 
(fTTjl is fulfilled, the left-hand side of (|S9^) is much larger 
than the right-hand side, and the electron-phonon cou- 
pling at small momentums can thus be safely neglected. 
As we mentioned earlier, for the Hubbard model near 
half-filling, the postulates (fTTj) breaks down, and the WB 
singularity is reachable; however, in this case the um- 
clapp effect becomes also important, and whether the 
WB singularity is robust against this additional effect 
remains an open question. Now we treat phonons as 3D. 
Plugging (H5"|) into (fB"T|) , we obtain 



Qtotal 



— E 

q 



2 , 2 ( ag 2 \ 4 



In- 



7 2 i (m_ 

It ~ \2n e 



(60) 

For small qs the negative retarded part in (f6"0"| is of the 
order q\ with a logarithmic correction, which is sub- 
dominant to q x . This implies that the WB singularity 
does not exist in 1D3D. This makes sense since the bulk 
phonon freedom suppresses phonon fluctuations on the 
chain and thus stabilizes the system against the attrac- 
tive electron-phonon coupling. Therefore, in a 1D3D sys- 
tem the electron-phonon coupling at small momentums is 
not important no matter the filling. A similar discussion 
for 1D2D leads to the same conclusion. 



V. A 3D ARRAY OF PARALLEL ID 
CONDUCTING CHAINS COUPLED TO 3D 
PHONONS 

In the previous two sections we have discussed two spe- 
cial systems: 1D1D and 1D3D. In this section we will 



8 



study a more realistic system: a 3D array of parallel ID 
conducting chains coupled to 3D phonons. 

In 1D1D we derived an exact relation between phonon 
and electron correlations, which is given in (|23|) . Per- 
forming an almost identical calculation we also find a 
similar result in 3D3D: 

( P 2kAr f )P2k F (0))+G ph (q), (61) 

where G^(q) = [g 2 +q 2 T + {gi/2n e ) 2 q 2 x + C ± ql}/iT 2 . We 
have replaced the sum over chain index j by the sum 
over rj_ with r± — nix (remember that the distance be- 
tween nearest-neighbor chains is tt after the dimensional 
rescaling given in (|11[) ). n an integer. We will adopt this 
change hereafter. 

For a further understanding of this problem, we follow 
the logic used in 1D1D and treat this 3D3D problem 
based on the same two limits. Let us start with the limit 
of weak coupling and/or high phonon frequency. In this 
limit we focus on the phonon-mediated effective electron 
action: 

S eff /h = -j± - 2g\ J <frj_d 2 r' ± dxdTdT' 

r±,r' ± 

{cos [2<f>(x, t, f±) - 24>(x, t' , f' x )] x 
G ph (T-T',f ± -r' ± )} , (62) 

where the phonon propagator is given by 



G p h(T, f±) 



d q±dq T 



q 2 + C ± q 2 ± + g{ 

(63) 



We have neglected the phonon dispersions along the 
chains, which only lead to insignificant correction. The 
phonon dispersions perpendicular to the chains, although 
typically very small, are important since they are respon- 
sible for the existence of long-range-order CDW. This will 
be seen later. Now we apply a self-consistent Gaussian 
variational method Q to this model . We approximate 
the action (f6"2"| as a quadratic one 



(64) 



Then we try to optimize the propagator G e} ^(q) so as to 
minimize the variational free energy 



F v — F + (S - So)c 



(65) 



where both the free energy F and the average ()q are 
calculated using the quadratic action Sq. Performing the 



procedure we obtain a self-consistent equation 

2 



G 



^2 f dT e" 4 ( <02>o " <0(o ^' 3)0(o ^^ )>o ) 



xG p/i (t, FjJ 



1 _ e i{q±-r±+tjTT) 



(66) 



We notice that the right-hand side of the above equation 
vanishes as q — > 0. Furthermore, if we naively expand 
the cosine in the action ((62)) to the quadratic order of 
4>, we will get corrections to |Vj_i?!>| 2 and (d T (f>) 2 and no 
corrections to (d x 4>) 2 . Therefore, we assume the following 
trial solution 

1 



G- el \q) = 



ir 3 K 



(C + Rgi 



Mql) 



(67) 



Hereafter we define K' as the renormalized K, given by 



K' = k/Vr. 



(68) 



The self-consistent equation (p6"|) can be solved analyt- 
ically in the limit M < 1, {K/K') 2 - 1 < 1. In this 
limit the system is highly anisotropic. Electron correla- 
tions along the chains decay as a power law at length 
scales shorter than 1 / \J~M and tend towards a non-zero 
constant at longer length scales, while electron correla- 
tions along _L directions decay very rapidly to the same 
non-zero constant. Specifically, they are given as 



-(<</» 2 )o-<0(O,O : O)0(O,T,r i )) o ) 
7--TT , f ± = 0,T < 1/VM 

Mt 1 , rl/0 or t> 1/VM 



(69) 



which reduce to the correlations of TLL in the limit M — > 
0. Plugging the formula (|69|) into the equation ([66]) we 
obtain two coupled self-consistent equations 



K 

K< 



-1 = 



i/Vm 



+32irKgl 



dT 



-2K'2 



1/VM 



G pH (t,0) 



1/VM 



M 



+1&-kKM k ' g 2 g^ . 
\QttKM k 'C±g%gi' 



where 



G ph (T,0) 



dT M K T 2 G ph (T,0) 

(70) 
(71) 

(72) 



-+al 



In the special case C± = 0, M vanishes, and the sec- 
ond and the third terms on the right-hand side of ([70]) 
thus vanish, and then (fTD)) reduces to the result of 1D1D. 
This is expected since the 3D3D model with C± = cor- 
responds to a simple sum of 1D1D models. For C± being 
nonzero, M can be easily calculated from (|7T|) as 



M = 



0. 



K' > 1 
K' < 1 



(73) 
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After a more involved, but essentially straightforward 
calculation, we get, to the lowest order of (K/K 1 ) 2 — 1, 



K 



2 2K-A 

929 1 



(74) 



for Cj_ <C <?? and 



K 
~K' 



M/C±, 



9hl K - 



M « g 2 1( , 
M » gf b) 



for C± ^> g 2 . The solution of M implies a quantum phase 
transition at about K = 1 for arbitrary weak electron- 
phonon coupling. This critical point can also be obtained 
by calculating the average of the retarded part of the ac- 
tion for a single chain with respected to the free electron 
action. This average scales as L 2 ~ 2K and thus also pre- 
dicts the critical point K = 1. For K < 1, electrons on 
different chains are effectively correlated, which leads to 
long-range-order CDW. The order parameter is the 2kp- 
wave component of the electron density, whose average 
is non-zero and scales as 



(P 



M 



2k r. 



2na 



(76) 



According to the general result given in (fBTj) . the 2k p- 
wave component of the phonon displacement field also 
becomes long-range-ordered with its expectation value 
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(77) 



For K > 1, electrons on different chains are effectively 
not correlated, and each chain behaves as an isolated 
single ID conducting chain coupled to 3D phonons, which 
has been discussed in section IIV1 For self-consistency, 
the solutions {721), d} and dZSJ) have to satisfy M < 1, 
(K/K') 2 — 1 <C 1, which leads to a restriction on the 
parameters 




C± < g 2 , 
C± » g\ ■ 



(78) 



This condition precludes the application of the Gaussian 
variational result in the region above OA in figure [5J 

For stronger coupling and/or lower phonon frequency, 
the result we get from the Gaussian variational method 
is not valid since the condition (|78p is violated. In this 
case we have to use the fixed-amplitude approximation, 
which has been discussed in section UTTl We assume that 
the important fluctuations of ip are its phase fluctuations, 
and the amplitude of ip can be approximated as fixed 
at its saturated value ip$. After this approximation we 



obtain a simplified model 

SfSjti = ^ J dxd 2 r ± dT^ [(d^) 2 + (d x ^) 2 } + 

(d T 0) 2 +g 2 (^j (d x 0) 2 +C ± \V ± 6\ 

xW>o}-^§^ J dxd 2 r ± dr cos (2<jy- 9) , 

(79) 

where we have transformed the action into its continuum 
version, and the constant pieces have been thrown away. 
We will calculate the fluctuations of </> and 6 based on 
this truncated model. 

The locally saturated value of ip is also estimated by 
using the mean-field theory, which approximates ip as 
uniform both in space and time. Since there are no bare 
electron correlations between different chains, the mean- 
field 3D3D model is a simple sum of the mean-field 1D1D 
models. Therefore, the mean-field calculation is the same 
as that in 1D1D, and the result is thus also the same; that 
is, ijiQ is nonzero and given by the formula (|29[) for K < 1, 
and zero for K > 1. However, unlike in ID ID, ip is now 
truly long-range-ordered for K < 1. In the following we 
will include the phase fluctuations and calculate the non- 
zero expectation value of tjj. 

Now we discuss the truncated model (|79]l . It can be 
shown that the cosine term in the model is always rele- 
vant no matter the value of ^q, as long as K is less than 2. 
The proof is given in appendix [Cl Therefore, for K < 1 
the fluctuations of (f> and 9 are always bound together, 
and both are characterized by the same action 

S[E)/h = J dxd 2 r ± dr [R(d T Z) 2 + B(8 X ~) 2 + 



MIV.SI- 



(80) 



where ^, can be either cf> or 9/2, parameters R, B and M 
are given by 



R 



(1 + AnK^l) . 



B = 1 + AnKgl 



2n f 



<Ao 2 



71/ 



(81) 
(82) 
(83) 



Note that the correction to B is negligible compared to 
the correction to R, since g\ is much less than 1. Plugging 
the expression of i/jq into (|8Tj) and (|83|) and writing R in 
terms of K' using (f6"8"]) , we get 



K 
~K 



-l^iglg 2 *-*) — 

M~C X (5 2 2 5^" 4 )^, 



(84) 
(85) 



which agrees with the result in reference [8J]. This result 
is consistent with the one obtained by the Gaussian vari- 
ational method in that both predict long-range-ordered 
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CDW and hence long-range-ordered ip for K < 1. The 
expectation values of p2k F an d V' are non-zero and given 
by the formulae (|76| and ([77]) . respectively. Clearly this 
fixed-amplitude approximation can not give result for at- 
tractive electron-electron interactions since the method 
relies on the condensation of the phonon field ip, which is 
only possible for repulsive electron-electron interactions 
according to the mean-field theory. Based on the result of 
weak coupling and/or high phonon frequency, we suspect 
that for strong coupling and/or low phonon frequency, 
there also exists a critical K above which both electron 
and phonon correlations become quasi long-ranged; how- 
ever, this critical K could be significantly bigger than 
1. 

Now let us discuss the validity of the fixed-amplitude 
approximation. Unlike in 1D1D, based on the above dis- 
cussions there is no self-consistent condition cause by 
phase fluctuations in 3D3D. However, the one caused 
by amplitude fluctuations still exists. After a similar 
calculation as we did in 1D1D, we find that the fixed- 
amplitude approximation is self-consistent only if 



fr » sT 
ff » £ 
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C±_ « g\ 
C± » g\ 



(86) 



which prohibits the application of the fixed-amplitude 
approximation in the region below OB in figure [2j 

According to the self-consistent conditions ([75)1 and 
(I86p . there is an intermediate region in the parameter 
space of gi and g2 in which neither of the approaches 
applies. This region is illustrated in figure [2] as the area 
between the loci OB and OA. In this region, we can not 
obtain a quantitatively trustable result; however, we do 
expect long-ranged electron and phonon correlations for 
repulsive electron-electron interactions. We also want to 
point out that for C± ~ gf, both approaches give the 
same result: 



M ~ I — 



fJi, 



(87) 
(88) 



This strongly suggests that in this special case, this quan- 
titative result should also hold in the intermediate region. 

So far we have focused on the electron-phonon cou- 
pling at large momentums. In section IIVI we discussed 
the effect of the electron-phonon coupling at small mo- 
mentums in 1D1D and 1D3D. Repeating the (virtually 
same) calculation, we find that the WB singularity also 
exists in 3D3D and is given by the formula ([55]) . Again, 
for a general filling, according to the postulates (fT"7|) the 
system is far away from the WB singularity; therefore, 
the electron-phonon coupling at small momentums is not 
important. 



VI. CONCLUSION 

In this paper we have studied electron and phonon cor- 
relations in systems of one-dimensional spinless electrons 
coupled to phonons at low temperatures. The focus is on 
the effect of backward electron-phonon scattering. We 
have been able to obtain quantitative results in the lim- 
its of weak coupling and/or high phonon frequency and 
strong coupling and/or low phonon frequency, which then 
leads to a qualitative understanding in the intermediate 
region as well. 

It has to be mentioned that our work for a three- 
dimensional array of one-dimensional conducting chains 
coupled to three-dimensional phonons does not take into 
account the effect of the interchain hopping. It is well 
known that the interchain hopping can break down the 
TLL behavior which we have assumed for electrons in 
our calculations. On the other hand, it is not clear if 
the effect of the interchain hopping will be suppressed 
in the presence of electron-phonon coupling. To under- 
stand these questions, a work which includes both effects 
of electron-phonon coupling and the interchain hopping 
will be interesting. 
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APPENDIX A: WEAK COUPLING AND/OR 
HIGH PHONON FREQUENCY IN 1D1D 

Since the total action is Gaussian in the phonon field tf>, 
we integrate out the phonon freedom and get an effective 
action only involving the electron field (f> 



SlP/h = 



^— J dxalT [{d T ^f + {d x 4>) 2 ] 
^ [ dxdrdr' GU(r - r') cos [2(f>(x, r) 

gi J 



2cj>{x,r')] . 



(Al) 



where the phonon correlation function is given by 

G ph {r)^e-^. (A2) 

We have neglect the phonon dispersions (d x <j)) 2 since 
their effect is not important. When we discuss an ar- 
ray of paralleled ID chains coupled to 3D phonons, we 
keep the phonon dispersions perpendicular to the wires 
since they can lead to long-range-order CDW. 

First we apply a Gaussian variational calculation to 
the phonon mediated effective action (|A1|) . We assume 
a Gaussian variational electron action 



S = J2 G el 1 (^*(m<l) 



(A3) 
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and try to optimize it to minimize the variational free 
energy 



F var — Fq + (S el — Sq)o j 



(A4) 



where Fo and ()o are calculated using action (|A3|) . After 
a straightforward calculation, we obtain 



<tf(fc,«r) = + fr) + —l dr e"*M x 

e -4[G e ,(o)- Gc! (r)] [! _ cos(g r r)] , (A5) 

where 

G el (r) = dq x dq T G el (q)e^ T . (A6) 

Since the right-hand side of the above equation vanishes 
as q — > 0, and the integral is independent of q x , we assume 
the following trial solution 



G el 1 (q x ,q T ) = ^J=q2 + Vtt£ 



(A7) 



with A' = K/yR. With the electron propagator writ- 
ten in this way, we ensure that the coefficient of q 2 is 
not renormalized but fixed at 1/ (wK ) . Hereafter we de- 
note A' and if as the renormalized and initial value re- 
spectively. Then the electron correlation function can be 
calculated as 



G el (0) ~ G el (r) = ^-ln 
Plugging this result into Eq. 



K' /K 2 ar + aK 



l 2 



aK' 2 
5l), we obtain 



(A8) 



1 

8<7 2 2 



A 2 



K 



— - 1 )lr 



dr e~ 9l|r| 

9l J-oo 

x [1 — cos(g T r)] . 



A 2 err + aK 



l 2 



aK 



l 2 



-2K' 



(A9) 



For q T — > 0, this leads to 



1 



A 



ttA V A 



'/ 2 
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+ 00 



dr T 2 e- 9llTl 



la 



aK t 
aK 12 



-2K' 



(A10) 



It is interesting to point out that in the absence of e _9l ' T ' 
in the integrant on the right-hand side of (|A10|) . the 
integral would converge for K' > 3/2 and diverge for 
A' < 3/2, and there would be a quantum phase transi- 
tion at the critical point K' = 3/2. The two phases are 
characterized by A' — > for A' < 3/2 and finite A' for 
A' > 3/2 respectively. With the presence of e~ 9l ' T ', the 
integral in (|A10p converges for any A'. Thus, there is no 



phase transition. Assuming A' is close to A, calculating 
K' to the lowest-order correction A' — A we get 



A 
A 7 



g- 2 \ n 2K-3 



- 1 



— In <7i 

91 yl 

Si 



A < 
A = 
A > 



(All) 



Although there is no phase transition, there exists a sharp 
crossover at A = K s for A < 3/2 with 



A' 



In (2a 
2 In 5 i 



(A12) 



which is obtained by doing (A/ A') —1 = 1. The effect 
of electron-phonon coupling on electrons is significant for 
K < K s and negligible for A > A s . This can be easily 
seen if we examine the renormalized plasma velocity v' — 
v j y/~R where R is given by 



R 



1 ||) gf<- 



(A13) 



Since both g\jg\ and g\ are typically much less than 1, 
for A > K s R is always close to 1, and v' gets hardly 
renormalized. For A — > A+, A suddenly becomes larger 
than 1, (Strictly speaking, formula (|A13[) doesn't apply 
anymore for A < A s , but it qualitatively tells the trends 
nonetheless.) and v' is thus reduced appreciably. In prin- 
ciple this result can be tested by experiments since K s 
can be less than 1 according to (|A12|) . 

Now we apply the standard momentum shell renormal- 
ization group calculation (RG) to this model. We sepa- 
rate the electron field into fast and slow varying com- 
ponents 4> = < / )> + < i where <j) > has support in the 
momentum shell ae di /a < q x < a/a, — oo < q T < oo, 
integrate out the fast varying component > and rescale 
the length with x — x'e M , t — r'e de so as to restore 
the original ultraviolet cutoff. The integration over <p > is 
achieved perturbatively by expanding the partition func- 
tion to the first order around the Gaussian fixed point. 
Performing the above procedure and using the expression 
(|A7|) for the propagator, we obtain the following RG flow 
equations (up to one loop) 



d[K'(i)] 





d£ 


d 


'&(£) 

gi V /_ 




di 


d[gx{l)] 


di 



-C 1 K r3 ^ 2 (e)f(g 1 (£)). 
9i 



(3-2AT')— (£), 
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(A14) 

(A15) 
(A16) 



where C± is a constant of order 1, f(x) is a step function 
which is 1 for x < 1 and for x > 1. We use f(x) to 
cut the renormalization of K'{£) at £\ — —\n(gia/a), 
where the phonon propagator G p /j(t — r') becomes effec- 
tively S(t — t'). The correspondence between our nota- 
tion and that of reference Q is: Y p h = 91/92, 9i ~ c^/^ph, 
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exp(2£) = K' . It seems that the correct eigenvalue in 
the RG flow equation (3.14b) in reference U should be 
2 — 2exp(2£ i ). Also there are some minor differences be- 
tween our RG flow equations and those in reference [||. 
They are due to different RG procedures and should not 
affect the physics. The RG flow Eq. (|A15|) seems to im- 
ply a phase transition at K' — 3/2; however, since the 
renormalization of K'(£) is utterly cut at £ = £%, K' is 
always finite even for K' < 3/2, so there will be no phase 
transition at K' = 3/2. For gi9\ 2 < 1, the perturba- 
tive RG holds way down to I = l\ since g\{£\) I g\{£\) 
remains less than 1, and we can integrate the flow equa- 
tions approximately to calculate K' as 



(K' - K) 



9i Jo 



(A17) 



which essentially gives the Gaussian variational result 
(jAllj) . For g 2 g?~ 2 > 1, siWMW becomes larger than 
1 at I = In [gi(£)/g 2 {£)] /(3-2K), which is less than l x . 
In this case the perturbative RG breaks down before we 
can calculate K'; therefore, a different strategy needs to 
be used to attack the problem. 



APPENDIX B: CALCULATION OF ^ 

After the mean-field approximation we get a simplified 
model 

qlD ClD -, /■ i 

-252^0 J dxdr cos (2</>) , (Bl) 
Stf/h =\J dxdr g\^l . (B2) 
The partition function can be calculated as 
Z = J D[MDW\ e -*<^+^+0 

- J D[ih] e- s l°/ h J D[<j>] e -^+0 

= [d[iPo] e-& s $ +F «U*» , (B3) 



where F e i[ipo] is the free energy of electrons for action 
SIP 4 
Gordc 
for ip 



Slf + S 1 ^ at a fixed ip . S^f + is exactly a Sine- 
Gordon model. Formula (IB3I) defines an effective action 



S; f h f bPo] = S 1 p Z + F el [ip ] 



(B4) 



To calculate F e i[ipo], we use Gaussian variational method 
9] . For any action So , we have 

F el [^]/h = ^ - hi (e-*(Si D +^ D -So) )o . (B5) 



where the () G denotes the average usin g th e action So- 
Given the convexity of the exponential [1J] one has al- 
ways 



(B6) 



and thus 

Flf < F var = F Q + (S- S ) Q . (B7) 
The basic idea is to optimize So such that F var gets to 

el 



F^ D as close as possible. Assuming that So is Gaussian 



Solh = \Y,^-q)G-\q)m 



with 



G-\q) = 



irK 



(B8) 



(B9) 



F var can be calculated as 



F var = -E hG ® + ^E« 2G ®- 

q>0 q 



2g 2 ipo \ dxdr e 



(BIO) 



Minimizing F var with respect to G{q) by doing 
dF var /dG(q) = leads to a self-consistent equation for 
D 



D 2 = 8TrKg 2( j > o 



D 



1 + Vl + D 2 



K 



(Bll) 



Assuming that g 2 ipo is much less than 1 (i.e., the single 
electron gap is much less than the Fermi energy) , we get 



D = 2 {2irKg 2 -4>o) 



l/(2-lf) 



(B12) 



for K < 2 and D = for K > 2. Therefore, K = 2 
is a naive estimation of the critical value of K at which 
the KBT transition happens for the Sine-Gordon model 
(|BTJ). Plugging (|B12|) into the formula (|BT0|) , we get the 
approximated F^ D 



i 



2ttK 



dxdr 



-(2 -if) (2<xKg2l>o)*=x 

(B13) 



where we have thrown out terms which are independent 
of ipo- Thus, according to (|B4|) the effective action for ifjo 
is given by 

ieff 1 



dxdr 



2-kK , 



(2 -if) {2irKg 2 ^o)~ + 
(B14) 



Minimize it with respect to , we get ipo — for K > 1 
and 



V-'o = Si 



2(7riT) 



9i 



(B15) 



for A" < 1. 
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APPENDIX C: CALCULATING PHASE 
FLUCTUATIONS FROM THE 
FIXED-AMPLITUDE MODELS 



We start with the truncated model (J26J). For conve- 
nience, we rewrite it as 

-2g 2 ipo J dxdr cos (20 - 9) , (CI 



where 



irK 



2 I x \ 2,2 

9i (j-J 1 X + 1, 



i% ■ (C2) 



Then we define two new variables, so that the action in 
terms of these two variables are decoupled. That is, the 
action can be expressed as 

S[T,n/h = ^G^(5)E( ? lS(- ? l + ^G^ l (9)x 
T(q)T(—q) — 252^0 / dxdr cosT, 



where 



T(q) = 2<f>(q) 



G^ 1 (^+4G, 1 ( (? ) 



G^iq) = , , 

G- 1 (< ? )+4G- 1 ( g ) 

G^ = GT 1 (q) + 4G e - 1 (g). 



(C3) 

(C4) 
(C5) 

(C6) 
(C7) 



Since T and £ are decoupled, using (|C4[) and (|C5|) the 
fluctuations of and <6 can be calculated as 



GT 1 ( 9 )+4G e - 1 (g1 



2 ' 



(C8) 



(9(q)0(-Q)) = 4G s (g) + 



G7 2 ( g -)(T( g -)T(-g)) 



Gl 1 (g)+4G e - 1 (< ? ) 



(C9) 



To calculate (T(g)T(— g)), we need to know whether the 
cosine term is relevant or irrelevant. Thus, we evaluate 



dxdr (cosT) = / dxdr e 2 , (CIO) 



where ()o is the average with respect to the harmonic 
part of the action (|C3p . Let us first calculate 



(T 2 ) n = 



1 



dq x dq T G r (q) 



dq x dq T 



0o 2 



(2tt)2 

1 

4ttK 
( Ue 

Plugging it into (|C10j) we obtain 

2-K- 

dxdr (cos T)o ~ L 



+ 



2K lni. 



(Cll) 



(C12) 



which implies a critical line in the ij} Q 2 -K plane defined 
by 



2is: = 4 . 



(C13) 



This critical line is illustrated in Fig. [T] The cosine term 
in the action (|C3[) is relevant below the critical line and 
irrelevant above it. In the two limiting cases the function 
(|C13|) predicts that ipfi — > oo for K — 2 and i/ifi ~ g 1 
for if = 0. For K < 1, it predicts 



4>o ~ 9i 



1/2 



(C14) 



For ipo < i/jq, the cosine term is irrelevant, and we get 
(T(q)T(-q)} = G r (q), which combined with (OH) and 
(IC9l leads to 



((l>m(-q)) = G^q) 
(0{q)9{-q)) = Ge{q). 



(C15) 
(C16) 



This implies that 9 and <f> are effectively decoupled. For 
ipo > V'oi C""( < ?)'^'( — 5)) i s massive, and we obtain 



(Hm-q)) = Gv(q), 
(9(q)9(-q)) = 4Gs(g) , 



(C17) 
(C18) 



which implies that the fluctuations of 9 and </> are bound 
together. 

Now we discuss the model (1751) . We follow the same 
routine as the above. Defining T and S using (|C4|) and 
(|C5j) . we get a decoupled action: 



S[T,E]/fi = ±J2 G iHoimX(-q) + lY,G^(q) 

252^0 



dxd r±dr cos T , 
(C19) 
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where and Gj}(q) are denned in ((C6j) and (|C7| . 

respectively. GT 1 (g) and Gg l (q) are now denned as 



we obtain 



TT 3 K 



(dl + 9?) , 



(C20) 



(C21) 



To check whether the cosine term is relevant or not, we 
evaluate the average of the cosine term for a single chain 
with respect to the harmonic part of the action. After 
we calculate 



(T 2 )o 



1 



(2tt)4 

1 



dq x d 2 q±dq T Gr(q) 



dq x d 2 q±dq T 



4n 3 K 



7 + 



2K\uL. 



(C22) 



dxdr (cos T)o ~ L 



2-K 



(C23) 



This implies that for K < 2, the cosine term in the action 
(|79[) is always relevant no matter the value of ipo . There- 
fore, for K < 1 and arbitrary ipo, the action (|T5|) predicts 
that the fluctuations of 4> and are bound together and 
given by (|C17p and (|C18p . respectively. 
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